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Abstract 

We review the physics potential for experiments with intense neutrino beams at the 
front-end of a muon storage ring, stressing the way the spin and flavour structure 
of neutrino interactions with matter can be used to shed light on the structure of 
the strong and electroweak interactions. Specifically, we discuss precision tests of 
the standard model, studies of polarized and unpolarized structure functions of the 
nucleon, and several new and exotic physics items. 
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1 Probing Matter with Intense Neutrino Beams 



Front-end physics at a neutrino factory is based on the realization that, be- 
cause of the flavour and spin structure of the coupling of neutrinos to weak 
currents, a neutrino beam is a unique probe of the structure of the standard 
model and of the structure of the nucleon. A neutrino beam thus has a greater 
physics potential than conventional electron or muon beams, provided the in- 
tensity of the beam is high enough. An accurate assessment of the physics po- 
tential of the experiments which could be performed with an intense neutrino 
beam cannot abstract from the fact that the time scale for the construction 
of a neutrino factory is of order of ten years: clearly, it is difficult to envisage 
what the 'standard model' will be ten years from now. Here, we will discuss 
the physics case for these experiments based on present-day knowledge. This 
is interesting not only because many of the measurements that we will discuss 
would only be possible at a neutrino factory (for instance, those related to 
polarized parton distributions), but also because comparison with what we 
already know will allow us to highlight the peculiar features of physics with 
neutrino beams. 

This brief review is largely based on a recent detailed quantitative study per- 
formed by a CERN working group [1]; quantitative estimates given here are 
taken from there unless otherwise stated. Previous studies on the physics po- 
tential of neutrino factories have been performed by working groups at Fer- 
milab [2] and Brookhaven [3]. General background on neutrino experiments 
is in Ref. [4]. Even though we will usually describe the energy and luminosity 
dependence of our results, we will assume the 'CERN reference scenario' [5]: 
specifically, a 50 GeV ji beam, with 10 20 muon decays per year along a 100 m 
straight section. 



2 Neutrino interactions with matter 



The most interesting neutrino-induced reactions in matter are neutrino-elec- 
tron elastic scattering, and neutrino-nucleon (deep)-inelastic scattering (DIS), 
i.e. neutrino-quark scattering. The former is a clean purely weak interaction 
process, while the latter allows one to probe strongly-interacting matter with 
weak currents. 



The cross section for elastic neutrino- or antineutrino-electron scattering is 

da 2G 2 F m e E v 



dy 



7T 



9i + 9R^-yy 



(1) 



1 



where y = E e /E u , < y < 1, so the total cross section is obtained re- 
placing (1 — y) 2 — > |. The couplings for neutral-current (NC) processes are 
9l = \{gv + Kg a), 9r = \{gv - Kg a), with \ u = -1 (A p = 1), and for 
charged-current (CC) processes g^ = 1 (<7l = 0), gn — = 1) for neutri- 
nos (antineutrinos). The total couplings are listed in the following table (the 
numerical values are computed with sin 2 9w = 0.23): 
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The total cross-section is tiny, of order a ~ 10~ 3 x (g 2 L + \g 2 R ) pb for 50 GeV 
neutrinos. Yet at a neutrino factory with a 20 ton liquid argon TPC or a fully 
active 2 ton liquid methane target one expects integrated luminosities of order 
of 8.6 • 10 10 pb _1 , leading to rates of order of ~ 10 7 events per year with a fi + 
beam (half with a beam). 

The neutrino-nucleon DIS cross section is by a factor ~ m p /m e larger: for 
charged-current interactions, up to corrections suppressed by powers ofm 2 /Q 2 
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where A are the lepton and proton helicities (assuming longitudinal proton 
polarization), and the kinematic variables are y — ^ (lepton fractional energy 

loss), x = (Bjorken x). The neutral-current cross-section is found from 
eq. (2) by letting mw — > rnz and multiplying by an overall factor [\{gv — 

Kg A )] 2 . ' 

At a neutrino factory, structure functions could be measured for 0.01 < 
x < 1 and 1 < Q 2 < 100 GeV 2 (note the kinematic limit s = ^ < 
2m p E fl = 100 GeV 2 ). Because y = Q 2 / \2xm p E v ) , at fixed x and Q 2 , y only 
varies with the neutrino energy. At a neutrino factory with a broad-band 
beam, neutrinos of various energies (measurable on an event-by-event ba- 
sis) are available. It is then possible to disentangle all structure functions Fj 
(unpolarized) or g^ (polarized) by fitting the y dependence of the data for 
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fixed x and Q 2 . One expects statistical errors of order of 1% or better on all 
three structure functions for x > 0.1 when Q 2 < 15 GeV 2 , and for x > 0.3 when 
Q 2 < 80 GeV 2 , and of order 10% otherwise. 

Information on the structure of the nucleon target is encoded in the structure 
functions, whose leading parton content in terms of the unpolarized and po- 
larized quark distribution for the i-th flavor ft = qp + qj^ and Aft = qj^ — qj^ 
is summarized in the following table, where for comparison we also give the 
standard result for charged-lepton scattering via virtual photon exchange: 

NC F? = \ Ei e 2 (ft + ft) 9l = \E t ef (Aft + Aft) 

NC F? = ±Zi(9v + 9A)i(Qi + Qi) 9i = \ Eifav + 9% (Mi + Aft) 

NC if = 2 Ei(gv9A)i (Qi + ft) 9i = ~ Ei(gv9A)i (Aft + Aft) 

CC F™ + =u+d+s+c g™ + = Au + Ad + As + Ac 

CC -F^ + /2 = u-d~s + c gf + = Au - Ad- As + Ac 

F 2 = 2xF x g A = 2xg 5 



Here ei are the electric charges and (gv)i, (9a)% are the weak charges of the 
i-th quark flavour. If W + — > (incoming P beam), then u <-> rf, c <-> s. 
The structure functions F 3 , #4 and #5 are parity- violating, and therefore not 
accessible in virtual photon scattering. Of course, beyond leading order in the 
strong coupling each quark or antiquark flavor's contribution receives 0(a s ) 
corrections proportional to itself and to all other quark, antiquark and gluon 
distributions. This last correction is flavor-blind, and thus decouples from the 
parity-violating structure functions F 3 , g± and g 5 . 



3 Tests of the Standard Model 



3.1 The weak mixing angle 



The weak mixing angle can be determined from the leptonic weak couplings 
(?£, and gn, which can be cleanly extracted from the measurement of the elas- 
tic cross-section eq. (1). The main background for this process is quasielastic 
neutrino-nucleon scattering. Because the transverse momentum of the outgo- 
ing electron is p t ~ ^/m e E v , but p t ~ ^m p E v for scattering off nucleons, the 
background can be removed with a p t cut. The current best determination 
from this process is 

sin 2 6 W = 0.2324 ± 0.0058 (stat) ± 0.0059 (syst) (3) 
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from the CHARMII experiment [6]. 
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Fig. 1. (a) Statistical uncertainty (in units 10~ 4 ) in the extraction of sin 2 9w from ve 
scattering as a function of the minimum electron energy, (b) Impact of luminosity 
measurement at the level of 10~ 3 , 5 10~ 4 , 1 • 10~ 4 on the same sin 2 6\y sensitivities. 

At a neutrino factory, events without a muon in the final state originating 
from or u e (in a /i~ beam) cannot be distinguished and must be considered 
together. With the luminosity discussed above, one then gets the statistical 
accuracy shown in Fig. la. The main systematic uncertainty is the determi- 
nation of the incoming neutrino flux (see Fig. lb), which can be done by 
normalizing to the inverse muon decay, namely the process u e e~ — > jj, 



or 



An independent determination of sin 2 9\y can be obtained from DIS. This has 
the advantage of larger statistics, but the disadvantage of systematic uncer- 
tainties related to knowledge of the nucleon structure. The value of sin 2 9w is 
extracted from the quark weak couplings: in ratios of NC to CC DIS cross- 
sections eq. (2) the lepton couplings cancel. Also, such ratios are less sensitive 
to details of nucleon structure, in that the leading-order dependence on parton 
distributions can be made to cancel. In particular, the Paschos-Wolfenstein [7] 
ratio 



— sin 2 9w 



(4) 



has been used to obtain a determination [81 



sin 2 w (O$) = 0.2255 ± 0.0018(stat) ± O.OOlO(syst). 



(5) 



At a neutrino factory, one can only measure combinations where and v e 
do not have to be disentangled on an event-by-event basis, such as W = 



°»cM+*Ncfr) R »+ = ^W+^W Qr p = °Nc(p-)-VNc(^) Qr Unear com _ 

bination thereof. These combinations give a good handle on sin 2 9 W , but are 
more dependent on nucleon structure: in fact, it is convenient to construct 
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an optimal combination, which maximizes the dependence on sin 2 9w while 
minimizing the uncertainty due to parton distributions. The uncertainties on 
sin 2 9w in units of 10~ 4 for various combinations are given in the following 
table: 



observable 


stat. error 


PDF 
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~ 12 


R» + 
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~ 15 


R»~ - 0.8iK 


2.2 


~ 2 


P 


4.9 


r~j 4 



In summary, at a neutrino factory the weak mixing angle could be mea- 
sured to an accuracy of about 10~ 4 , in two different ways. This is com- 
parable to the best presently available determination of sin 2 ^, namely [9] 
sin 2 6w = 0.23098 ±0.00026 from all available asymmetry measurements. Fur- 
thermore, the two measurements from neutrino elastic scattering and DIS 
undergo different radiative corrections from each other and from the asym- 
metry measurements. This would therefore be a very competitive test of the 
standard model, or its violation. 



3.2 The strong coupling 



The strong coupling can be extracted from DIS data either by considering 
combinations which do not depend on parton structure (sum rules), or by 
performing global fits where both a s and all parton distributions are simulta- 
neously determined. In both cases neutrino beams are superior to conventional 
charged lepton beams because of the availability of more independent combi- 
nations of individual parton distributions. 

Specifically, if both proton and neutron targets are available, one may con- 
struct two combinations which only depend on the total number of up plus 
down valence quarks (Gross-Lewellyn Smith sum rule [10]) or up minus down 
valence quarks (unpolarized Bjorken sum rule [11]), up to an a s -dependent 
factor which is currently known to 0(af): 
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where the coefficients of the power corrections /igls an d /^Bju are unknown and 
must be fitted. The disadvantage of this way of determining a s is the need to 
extrapolate over the full range < x < 1 the data which are only available 
in a limited range of x. The accessible x range is larger at lower Q 2 , where 
however power suppressed corrections are larger. At present, a determination 
of a s (Mz) with an error Aa s (Mz) = looi2 has been obtained from the GLS 
sum rule with the CCFR neutrino beam [12], while the BjU integral has never 
been measured because of the impossibility of disentangling F 1 and F 2 from 
present-day neutrino DIS data. At the neutrino factory, one could reach an 
uncertainty Aa s (M z ) = 0.0035, which is completely dominated by limited 
kinematic coverage in either Q 2 (power corrections) or x (small x extrapola- 
tion) and could only be improved if a higher energy beam were available. 

A very competitive determination could be obtained by performing a global 
fit to structure functions: interestingly, because more structure functions are 
available, at a neutrino factory the full set of parton distributions could be 
determined in a single experiment (see Sect. 4 below). One could then achieve 
a statistical accuracy on a s of order Aa s (M z ) = 0.0003, to be compared to the 
statistical accuracy Aa s (M z ) ~ 0.0017 which can be obtained from present- 
day global fits [13] which include both v and charged-lepton DIS data from 
various experiments (the error from current neutrino experiments alone is of 
order Aa s (Mz) ~ 0.005 [14]). The statistical precision of this determination of 
a s value is by one order of magnitude better than the extant global error [15] 
Aa s (M z ) ~ 0.003. The accuracy in the determination of a s at a neutrino 
factory would thus be entirely dominated by theoretical uncertainties, and it 
could be the most accurate determination once next-to-next-to-leading order 
corrections to perturbative evolution [16] are known. 

More studies of strong interaction physics could take advantage of the fact that 
charm is copiously produces in charged-current events and easily detected. 
This could be exploited not only in order to further refine our knowledge 
of parton distributions through the study of specific semi-inclusive channels, 
but also to study QCD corrections to charm production near threshold, and 
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finally to measure accurately absolute branching ratios and decay constants 
of individual charmed mesons, such as A c or D s . 



4 The Structure of the Nucleon 

Current information on the parton distributions of the nucleon [17] comes 
mostly from DIS data, the bulk of which are produced with charged lepton 
beams and thus are essentially NC scattering. The problem is then that (see 
Sect. 2) only one combination of qi + q\ distributions is accessible. This means 
that different quark flavors can only be disentangled using isospin, if proton 
and neutron targets are available, plus in principle by exploiting perturbative 
evolution (i.e. subleading corrections). Hence, using NC DIS data it is hard 
to determine strangeness, and impossible to measure the C-odd combination 
q-q. 

This is to be contrasted with the situation in CC scattering [18], where if 
both proton and neutron targets and v and v beams are available, then one 
can form eight linear combinations of the two independent structure functions 
(F 1 and F 3 unpolarized, gi and g 5 polarized), six of which are independent 
(NC data do not give any extra independent information). It is easy to con- 
struct leading-order combinations of parton distributions which, below charm 
threshold, determine all six light flavours and anti-flavours independently. The 
charm distribution can then be determined either comparing data below and 
above charm threshold, or tagging charm events, which have a distinct dimuon 
signature [4] . In practice, of course, parton distributions will be determined by 
fitting the full next-to-leading order expression of structure functions, however 
the fact that individual partons can already be disentangled at a leading-order 
guarantees the accuracy of the NLO determination. 

4-1 Unpolarized DIS: the flavor content of the nucleon 

Unpolarized parton distributions are a necessary ingredient in the computation 
of any collider process. However, only the up, down and gluon distributions 
can be determined in a reasonably accurate way from DIS data. Some infor- 
mation on strangeness can be extracted [19] from neutrino data, while some 
less-inclusive observables (such as W production, or Drell-Yan) provide some 
constraints on the relative size of the q and q distributions, but the results are 
at best semi-quantitative (see Fig. 2a). 

A quantitative estimate of the accuracy at a neutrino factory can be ob- 
tained by generating pseudo-data for structure functions with appropriate 
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Fig. 2. (a) Percentage error on parton distributions at a neutrino factory (solid) 
compared with present-day [13] errors (dashed), (b) Correlation coefficients between 
parton distributions determined at a neutrino factory. 

errors, and then producing a fit of parton distributions based on this data. 
In Fig. 2a the error estimates on individual partons obtained in this way are 
compared to the extant knowledge. Note that no current errors on strange 
and antiquark distributions are given, since the present results largely depend 
on theoretical prejudice. In Fig. 2b we further show that the point-by-point 
correlation of individual distributions is uniformly quite low, indicating that a 
model-independent flavour and antiflavour separation is possible to 10%-20% 
accuracy in most of the accessible kinematic range. A neutrino factory would 
essentially turn the determination of individual parton distributions into a 
precision quantitative exercise. 



4-2 Polarized DIS: the spin of the nucleon 



Polarized DIS has recently attracted considerable attention because the un- 
expected smallness of the proton's singlet axial charge suggests that the 
nucleon spin structure is considerably subtler than naive parton expectations 
might suggest [20]. In the naive parton model the singlet axial charge is the 
fraction of the nucleon spin which is carried by quarks; the Zweig rule would 
lead one to expect this to be around 60% but the experimental value is com- 
patible with zero. 

This state of affairs point to several possible scenarios for the nucleon spin 
structure. Beyond leading order the axial charge is given by 

a = AS - ^AG, (7) 
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where AS = J2i(Aqi + Aq~i) is the scale-invariant quark spin fraction, and AG 
is the gluon spin fraction. The latter, due to the axial anomaly, gives an effec- 
tively leading-order contribution to a 7 (it depends on scale as AG ~ ^-). 
So a first possibility is that AG is large enough that the quark spin AE is 
large even though the axial charge a is small ('anomaly' scenario [21]). A 
different option ('instanton' scenario) is that AE is itself small, because of 
a large contribution from sea quarks whose polarization is anticorrelated to 
that of valence quarks (possibly because of 'instanton' QCD vacuum configu- 
rations [22]). Yet another possibility ('skyrmion' scenario) is that AE is small 
because of a large contribution from 'valence' strange quarks [As] >> |As| 
(as suggested [23] in the Skyrme model). 

At present, the quark and gluon spin fractions can be extracted from NC DIS 
data [24]: AG(1, 1 GeV 2 ) = 0.8 ± 0.2, AE(1) = 0.38 ± 0.03, while the 'octet' 
combination can be determined using SU(3) from baryon /5-decay constants: 
ag = Au + Ad — 2As = 0.6 ± 30% (the large error comes from a conservative 
estimate of SU(3) violation). It is of course impossible to polarize the kind 
of targets which are required for present-day neutrino DIS experiments, so 
no neutrino data are available, and thus little information on strangeness and 
no information at all on the quark-antiquark separation is available in the 
polarized case. At a neutrino factory, significant rates could be achieved with 
small targets [18]: with a detector radius of 50 cm, 100 m length, the structure 
functions g±, g§ could be independently measured to an accuracy which is 
about one order of magnitude better than that with which gi is determined 
in present charged lepton DIS experiments. 

We can then assess the information that one might obtain on the nucleon 
spin structure by generating pseudodata within different representative scenar- 
ios [25]. Assuming an 'anomaly' scenario one would determine Ag = 0.9 ±0.1; 
AE = 0.39 ± 0.01; a s = 0.56 ± 0.01, while in an 'instanton' scenario Ag = 
0.2 ± 0.1; AE = 0.32 ± 0.01; a s = 0.57 ± 0.01. Furthermore, in an 'instanton' 
scenario [As - As] (1, 1 GeV 2 ) = -0.007 ± 0.007; while in a 'skyrmion" sce- 
nario one would observe [As - As] (1, 1 GeV 2 ) = -0.106±0.008. Note that the 
strange contribution, and hence the octet component, would be determined 
from the data directly, without having to use SU(3) (which would thus be 
tested in the process). Clearly, the distinct scenarios could be well separated 
from each other: in fact full flavor separation at the level of first moments 
would be possible. Hence, a resolution of the nucleon spin structure would be 
possible at a neutrino factory. 

The situation would be less favorable for the x-dependence of parton distri- 
butions, because of the difficulty of disentangling the potentially large gluon 
contribution. If the gluon were well determined from other experiments, how- 
ever, a full point-by-point determination of the polarized parton content of 
the nucleon would be possible, to an accuracy comparable to that of the un- 
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polarized case. 



5 Exotica 



On top of the standard processes discussed so far, the high luminosity at a 
neutrino factory allows the study of several rare or exotic processes, of which 
we give two examples. 

a) b) 

^ V 

W" ■. c r 




D s 4 D 



N N N N 

Fig. 3. Feynman diagram for hard exclusive D s production 

A first example is the extraction of off-diagonal (or skewed) parton distribu- 
tions (SPD) from hard exclusive meson production. The Feynman diagrams 
for this process are shown in Fig. 3. If the virtuality Q 2 of the W boson is 
large compared to the nucleon momentum transfer t and all masses, then the 
cross section for the process factorizes [26] as dXB d dQ 2 dt (W~ + N — > D~ + N) = 
H ® $£> ® F, where H is the cross-section for the underlying hard perturba- 
tive parton subprocess, $£> is the D~ fragmentation function, and F(x,t,Q 2 ) 
is an SPD, which interpolates between the usual nucleon form factor G(t) 
(which is related to the first :r-moment of F) and parton distribution F(x) 
(related to the t — > limit of F). The estimated cross section for this process is 
a = 2.2 x 10~ 5 pb, leading to 10 4 event/yr at the neutrino factory. Given good 
knowledge of the D s fragmentation, the SPD can be measured. Because of rel- 
atively low backgrounds and the distinct experimental signature, this would 
be a favourable way of measuring SPDs, which have never been determined 
so far. 

As another interesting 'exotic' process, consider neutrino-electron annihilation 
into hadrons through radiation of an intermediate virtual W [3]. One can then 
define the analogue of i?-ratio of e + e— annihilation: 

„ a{v e e~ — > hadrons) . . 

Ra = —r 1 - 8 



Ra(s) is the spectral function for annihilation into axial vector final states: at 
large y/s, Ra(s) can be computed in perturbative QCD from the underlying 
parton processes, while at low y/s ~ m n , Ra(s) can be determined using 
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PCAC from the relevant axial current matrix element. The threshold energy 
for the inverse muon decay process (z^t/i~ final state) is £7^ = 11 GeV, while 
the threshold for the lowest hadronic final state 77r~ is E u = 19 GeV, so at a 
neutrino factory one could perform some interesting tests of PCAC and chiral 
symmetry. 



6 New Physics 

All the standard physics items which we discussed so far can also be viewed as 
tests for new physics, or tools in searches for new physics. Specifically, different 
determinations of the weak mixing angle lead to tests of the electroweak sector 
through the comparison of the corresponding radiative corrections (Sect. 3.1). 
Determinations of the strong coupling (Sect. 3.2) and especially its running 
test the strongly interacting sector. Precision determinations of parton dis- 
tributions (Sect. 4) are a necessary input in searches for new physics: in fact 
some recent possible indications of new physics could also be explained by 
invoking lack of accurate knowledge of parton distributions [27]. 

However, high-intensity neutrino beams can also be used to design specific 
searches for new physics, of which we also give two examples. First, one 
can exploit the copious production of charmed mesons to study specific non- 
standard decay channels. An interesting possibility is the search for possible 
T-violation in the decays of the A c [3]. Indeed, in the semileptonic decay 
A+ — > l + uA one can construct an experimentally observable T-odd corre- 
lation Ct = (<7a • {pa x Pi)) from the spin and momentum of the daughter 
hyperon and the momentum of the final-state lepton . This correlation can- 
not be affected by either strong or electromagnetic final state interactions, so 
a nonzero value for it is a measure of T-violation. 

Another example is the search for lepton-flavour violating decays, such as 
ji — > eu T u.i. The corresponding four-Fermi couplings could be probed at a 
neutrino factory down to a strength about ~ 3 x 1CT 4 weaker than ordinary 
weak interactions. Of course stronger bounds already exist in the charged 
lepton sector, but the current limits for processes involving neutrinos are only 
of order ~ KT 1 - 10 -2 . 



7 Outlook 

The peculiar features of a neutrino beam is the availability of a probe which 
depends both on spin and flavor. However, the physics potential of such a 
beam can only be exploited given high enough intensity. Consequently, there 
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is a whole class of measurements which is only possible at a neutrino factory. 
An example is charged-current polarized deep-inelastic scattering, which can 
only be studied accurately with a neutrino beam and a reasonably sized target: 
with charged lepton beams it is hard to separate charged-current and neutral- 
current events, and with low intensities it is impossible to polarize the target. 
Which of these measurements might be the most interesting at a neutrino 
factory will largely depend on the development of high-energy physics in the 
next decade. In any case, a neutrino factory has the potential of developing a 
broad program of short-baseline physics, which could answer many outstand- 
ing questions which are not likely to be addressed at any other experimental 
facility. 
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